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Abstract 

We report resonant x-ray scattering data of the orbital and magnetic ordering at low tempera- 
tures at the Mn £2,3 edges in Lao.sSri.sMnO^ The orderings display complex energy features close 
to the Mn absorption edges. Systematic modeling with atomic multiplet crystal field calculations 
was used to extract meaningful information regarding the interplay of spin, orbital and Jahn- Teller 
order. These calculations provide a good general agreement with the observed energy dependence 
of the scattered intensity for a dominant orbital ordering of the d x 2_ z 2 / d y 2_ z 2 type. In addition, 
the origins of various spectral features are identified. The temperature dependence of the orbital 
and magnetic ordering was measured and displays a strong interplay between the magnetic and 
orbital order parameters. 

PACS numbers: 61.10.-i,71.30.+h,75.25.+z,75.47.Lx 



I. INTRODUCTION 



Orbital ordering, which involves correlations between the spatial distribution of the outer- 
most valence electrons, has long been considered a vital ingredient in the structural and phys- 
ical properties of strongly correlated electron systems such as transition metal oxides. The 
competition and cooperation between the charge, orbital, and spin degrees of freedom of the 
electrons manifest itself in unusual properties such as high temperature superconductivity, 
colossal magnetoresistance and magnetostructural transitions^. In particular, charge-orbital 
ordering in half-doped manganites has attracted much attention and controversy ^ 3 ' 4 ' 5 ' 6 . In 
Lao.5Sr L5 Mn04 the Mn sites are all crystallographically equivalent at room temperature 
with an average valency of +3.5. This material displays a phase transition at ~ 240 K 
below which it was believed that charge disproportionation of the Mn ions occurs^, cre- 
ating two inequivalent sites identified as Mn 3+ and Mn 4+ . This was originally predicted 
by GoodenouglA It should be noted that the picture of integral charge on the Mn ions 
has recently been challenged^, implying that the actual valence difference of the two Mn 
ions is much less then 1. The assumption that they are close to 3.5+, however, is based 
on the formal valence and does not take covalency into account which should cause a re- 
duction; according to Hartree-Fock calculations on the related compound Lao.sCao.sMnOj^ 1 - 
the charge on the manganese ion is even below 3+. In this paper, we adopt a multiplet 
picture which applies to integer charge only and we assume the value of 3+ as the closest 
to the real situation. The strong Hund's rule coupling and the cubic (Oh) component of 
the crystal field implies that the Mn 3+ 3d 4 site has three electrons in the level an d one 
electron in the twofold degenerate e g ^ level. The degeneracy of the e g ^ level can be lifted by 
cooperative Jahn- Teller distortions of the Mn06 octahedra reducing the symmetry to D^- 
Adding the tetragonal field, ti g \ level splits into d xy and doubly-degenerate d xz>yz orbitals, 
while e g i level splits into d 32 2_ r 2 and d x 2_ y 2. Therefore, depending on the value of the ap- 
plied tetragonal crystal field, Mn 3+ will have as a ground state filled t 29 T levels and the 
orbital lying along the elongation direction, either d 3z 2_ r 2 or d x 2_ y 2 orbital (see Fig^ and^3 
respectively), which will determine the type of the orbital ordering. Further cooling below 
Ttv ~ 120 K results in long range anti-ferromagnetic ordering of manganese ions into a CE 
type structure 1 ^. This was the basis of an alternative model put forward by Goodenough^ in 
which orbital order correlations lower the energy by favoring anti-ferromagnetic interactions. 



a) 





FIG. 1: The structure of Lao.sSri.sMnO^. a) and b) show two different types of orbital ordering, 
d 3x 2_ r 2/d 3y 2_ r 2 type in a) and d x 2_ z 2 / 'd y 2_ z 2 in b). c) shows the magnetic ordering with blue lines 
representing the ferromagnetic spin chains. The orbitals are drawn for the Mn 3+ ions, while the 
small yellow circles represent oxygen atoms and large blue circles represent the Mn 4+ ions. 



Until very recently, detailed investigations of the origin of the orbital order remained elusive 
because of the lack of a technique capable of direct observation. 

Resonant x-ray scattering (RXS) studies of Lao.5Sr L5 Mn0 4 were first attempted at the 
manganese K edge. The diffraction intensity at an orbital ordering reflection displayed a 
striking resonant enhancement near the absorption edge and a specific dependence on the 
azimuthal angled. However, the Mn K edge resonance involves virtual excitations from Is 
to Ap bands and is thus generally insensitive to orbital ordering of the 3d states. Theo- 
retical studies proposed that the observed sensitivity is largely due to Jahn- Teller distor- 
tions and 4p band structure effect a 14 ! 15 ! 16 ' 17 ' 18 ' 19 , rather than 3d — 4p Coulomb interactions^ 
and this conclusion was recently supported by experimental work on strained manganite 
films^i. Experiments performed at the L edges, on the other hand, directly probe the 3d 
states. Resonant soft x-ray diffraction at the Mn L edges was first reported by Wilkins 
et al^ and was followed by the first direct observations of orbital ordering using soft x- 
rays in Lao.sSri.sMnO^. Aided by earlier theoretical predictions^ they reported energy 
resonances at the Mn L 3 and L 2 edges showing that the orbital ordering was caused by a 
mixture of both cooperative Jahn- Teller distortions and direct Goodenough orbital correla- 
tions. Similar experimental results were later independently reported by another group^. 
More recently, the same technique has been used to study magnetic and orbital correlations 
in Pr . 6 Cao.4MnO^. 

In this paper, we report both experimental and theoretical results on resonant x-ray 
scattering from orbital and magnetic order at low temperatures in Lao.sSrx sMnO^. We 
describe the measured intensity of the specific peaks in the energy spectra for both the 
magnetic and orbital reflections. We also report calculations using an atomic multiplet 
picture incorporating the crystal field effects. We present the spectra obtained by fitting the 
crystal field parameters and discuss their agreement with the experimental spectra. Further, 
we indicate the dominant orbital ordering type. We also present analysis of the superlattice 
reflections in terms of the Jahn- Teller vs. orbital ordering origin of the spectral features. 
Finally, we have measured the temperature dependence of certain features in the energy 
resonant spectrum. These show a dramatic change at a temperature corresponding to the 
Neel temperature, indicative of a strong interplay between the magnetic and orbital order 
parameters. 

This paper is organized as follows: Section II of the paper describes the methods to obtain 



the spectra, experimental and theoretical. The results are presented in section III, where 
experimental spectra and theoretical fits are compared. There we also present the analysis 
of the Jahn- Teller contributions to the orbital ordering scattering. Section IV describes 
the temperature dependence of the spectra, while all our results are summarized in the 
conclusions. 

II. METHOD 

A. Experimental approach 

The experiments were performed at beamlines 5U1 at Daresbury and ID08 at the Eu- 
ropean Synchrotron Radiation Facility (ESRF). Single crystals of Lao.sSrx 5 Mn04 with di- 
mensions 10 x 3 x 3 mm 3 were grown at the University of Oxford using the floating zone 
method. They were cut with either [110] or [112] directions surface normal and polished 
with 0.25 \xm diamond paste to a flat shiny surface. 

The crystals were mounted in the ID08 5-circle diffractometer operating at a base pressure 
of 1 x 10~ 8 mbar, and equipped with a Si diode detector. The beamline produces ~ 100 
% linearly polarized X-rays. At ~650 eV an energy resolution of 165 meV was obtained. 
Sample cooling was achieved using a liquid He cryostat attached to the sample stage by 
copper braids resulting in a base temperature of 63 K. 

On the 5U1 beamline, crystals were mounted on a two circle diffractometer enclosed in 
a high vacuum chamber with a base pressure of 1 x 10~ 8 mbar. The incident beam has a 
resolution of 500 meV, and a beamsize of 1 x 1 mm 2 . The sample was mounted on a liquid 
nitrogen cooled copper block, with a small manual x adjustment for initial orientation. The 
base temperature achieved was 83 K. Temperature stability between 83 K and 300 K was 
achieved with a heater element, and control thermometer mounted next to the sample. 

At both beamlines the experimental procedure was identical. The incident energy was 
set to the manganese L 3 edge, and superlattice peaks located at the (7,4,0) and (4,-4,0) 
positions. Each reflection was measured on a separate sample. Energy scans of the reflections 
were performed at fixed wavevector. The integrated intensity was measured by longitudinal 
scans through the peak at fixed energy. 

The azimuthal dependence of the orbital order superlattice reflection was measured on the 



ID08 diffractometer by rotating the sample around the scattering vector, q. In the absence 
of an orientation matrix (UB) this was achieved by rotation of the diffractometer's <fi axis. 
The azimuthal angle \& is 90° when the [001] direction lies within the scattering plane. 

B. Theoretical model 

The calculations of the RXS spectra are based on atomic multiplet calculations in a crystal 
field. Cowan's atomic multiplet progra m 27 i 28 i 29 provides ab-initio (Hartree-Fock) values of 
the radial Coulomb (Slater) integrals F ' 2, • (d,d), F 2 (p,d), G 1,3 (p, d) (direct and exchange 
contributions) and the spin-orbit interactions ((2p) and C(3rf) for an isolated Mn 3+ ion. 
Their values are given in table H] As we are lowering the symmetry by the inclusion of the 
crystal field, we apply the Wigner-Eckart theorem to calculate matrix elements in a given 
point group* 3 ^ starting from the matrix elements evaluated in the spherical group. This is 
implemented in the "Racah" code, which also gives the values of the dipole transition matrix 
elements, necessary for the calculation of the spectra, as will be described in the following 
subsection. The fitting procedure starts with adjusting the crystal field type and strength 
by modifying the cubic (X 400 ) and tetragonal (X 220 ) crystal field parameters. Orbital and 
magnetic spectra calculated using the atomic multiplet code and including the crystal field 
effects are then compared to the experimental ones and the procedure is continued until the 
optimized set of crystal field parameters is found. In order to take into account the screening 
effects present in a real crystal with respect to the atomic picture, we scale down all the 
Slater integrals to 75 % of their atomic values. 

In the actual Lao.5Sri.5Mn04 compound, the spins are approximately aligned in the [110] 
direction 4 ^, displaying anti-ferromagnetic ordering along the c-axis and between the spin- 
chains, and ferromagnetic along the spin-chains, as shown in Fig. In order to simulate 
the effects of the spin ordering on an isolated ion, (superexchange and direct exchange 
interactions) we introduce a magnetic field acting on the spin of the atom. This field splits 
additionally the S = 2 quintet into S z = —2, — 1, . . . , 2 levels. Inclusion of the magnetic field 
favors the S z = —2 level as the ground state of the atom. The strength of the exchange 
interaction was set to 0.02 eV. 

Using the above procedure one readily obtains the eigenvalues for the ground (3d n ) and 
excited (2p 5 3d n+1 ) states separately. The next step towards the calculation of the scattering 



F 2 (d,d) F*(d,d) F 2 (p,d) G\p,d) G 3 (p,d) C(2p) C(3d) 
3d 4 11.415 7.148 - 0.046 
2p 5 3d 5 12.210 7.649 6.988 5.179 2.945 7.467 0.059 

TABLE I: The Hartree-Fock values for the ground and excited states of Mn 3+ given in eV. Coulomb 
and exchange integrals in the calculation have been scaled to 75 % of their atomic values. The 
p-shell spin-orbit parameter C(2p) has been increased by 9 % from the Hartree-Fock value to 
correspond to the experimental valued. 



or absorption spectrum is the evaluation of the dipole transition probability. Depending on 
the crystal symmetry and the type of scattering one is interested in, the angular dependence 
and the expectation value of the transition operators of the resonant amplitude will vary^ 2 .. 

The crystal field has the tetragonal (D^) symmetry. However, if the spin direction is 
assumed along [110], the site symmetry is lowered to that of the Cj point group. Following 
the approach described by Carra and Thole^ 2 -, we derived the expression which describes 
the atomic resonant scattering amplitude^ and computed the structure factor. In the case 
of the orbital scattering with the wave vector (|,|,0), the resonant scattering intensity is 
proportional to the following combination of the atomic scattering tensor components: 

I°° oc \ F° — F° 4- F c — F° I 2 CFl 

J res ^ I r 10:1-1 r 10:11 ^ r 11:10 ^1-1:101 ' \ L ) 



with F* m . lm , defined as: 



lm;lm' ^ E - E n + Kuj + iT /2' 1 ' 



where m and m' denote polarization states and are the electric dipole operators denned 
in spherical coordinates. |0) represents the ground state with energy Eq and \n) intermediate 
state with energy E n . The photon energy is Tiu and V stands for the broadening due to the 
core-hole lifetime. 

Similarly, for the magnetic scattering with the wave vector (\,— \,\), the scattering in- 
tensity can be expressed as: 

T u0 oc \ F° - F e I 2 (X\ 

-'res ^ rilill ^l-ljl-ll ■ \°) 

This can be written using the more usual notation^ in which i^ s ° oc \F^ — -F 1 c _ 1 | 2 . 

In these calculations we used T = 0.5 eV and in addition, the scattering intensity was 
convoluted with a Gaussian of width 0.1 eV to simulate the experimental (energy) resolution. 
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FIG. 2: Energy scan through the (2,3, 0) orbital order reflection at constant wavevector at 63 K 
(full black line) with the theoretical fit (dashed red line). The inset shows the temperature evolution 
of the orbital order parameter. 

III. ORBITAL AND MAGNETIC SCATTERING 

Figure |2] shows the energy dependence of the scattered intensity at a fixed wavevector of 
Qoo = {\i\$) through the Mn L 3 and L 2 edges. This wavevector corresponds to the long 
range orbital ordering of the Mn 3+ ions within the a — b plane. These data were collected 
at the beamline ID08 at the ESRF, Grenoble, France, on a single crystal cut with the [110] 
direction surface normal. On initial inspection, the spectrum is dominated by scattering at 
the L 3 edge where two distinct features are present. The scattering is weaker at the L 2 edge 
where three convoluted peaks are observed. Figure El shows the energy dependence collected 
in a similar fashion, but at a wavevector of (Jaf = (i? - \^\)- This reflection corresponds 
to long range anti-ferromagnetic ordering on the Mn 3+ sub-lattice. The observation of a 
c-axis component in the wavevector of the anti-ferromagnetic reflection (compared to the 
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FIG. 3: Energy scan through the (4,-4,2) magnetic order reflection at constant wavevector at 
63 K (full black line) with the theoretical fit (dashed red line). The inset shows the temperature 
evolution of the magnetic order parameter. 

orbital order wavevector) indicates that whilst the orbitals are ferro-ordered along the [001] 
direction, the magnetic moments are aligned anti-ferromagneticaliy. The magnetic ordering 
energy spectra is even more dominated by scattering at the L 3 edge, which is observed to 
contain up to five separate features. The L 2 features are much weaker and contain only 
two main peaks. These data show no relative shift in energy between magnetic and orbital 
spectra, at variance to what was reported by Thomas et alM: The insets of Figures 121 and El 
show the temperature evolution of the orbital and magnetic order parameters respectively. 
The orbital order parameter was found to decrease in a continuous fashion with increasing 
temperature in the vicinity of Tqo = 230 K. The anti-ferromagnetic reflection shows the 
same behavior, albeit, with the lower transition temperature of T/v = 120 K. 

Figure |U displays the integrated intensity of the orbital order superlattice reflection as 
a function of the azimuthal angle The solid line shows a simulation of the azimuthal 
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FIG. 4: Dependence of the integrated intensity as a function of azimuthal angle ^ for the (|,|,0) 
superlattice reflection (circles). The solid line is a simulation of the azimuthal dependence for the 
proposed orbital order. 



dependence given by the following equation 

1(9, V>) oc cos 2 6 B cos 2 * • I r °°, (4) 

where 9b is the Bragg angle of the reflection and xjj is the azimuthal angle. This dependence 
is characteristic of the El resonant scattering process. It is worth mentioning that this 
azimuthal dependence, which shows an extinction of the scattering when \l/ « 90° is not 
consistent with scattering from the magnetic structure with moments aligned along [210]. 
Therefore we conclude that the scattering observed at the (|,|,0) position arises solely from 
the orbital ordering. 

In Fig. |21 we also present the best fit we have obtained for the orbital scattering by 
fitting simultaneously to the both data sets. The fitting procedure yielded the following 
optimized crystal field parameters: X 400 = 5.6 eV and X 220 = 3.75 eV, or 10D q = 1.70 eV 



and D s = —0.45 eV . This corresponds to orbital ordering of the d x 2_ z 2 / 'd y 2_ z 2 type, as 
illustrated in Fig ^p. This is in contrast to the conventional model of orbital ordering of 
the d,3 X 2_ r 2/d 3y 2_ r 2 typ o 14 i 24 i 25 , but in agreement with the recent findings of Huang et al.—. 
Magnetic scattering corresponding to the same crystal field parameters is shown in Fig. 
01 All the features of the L3 and L2 structures of the orbital scattering fit display good 
agreement with the experimental data although the L 3 /L 2 branching ratio is somewhat 
overestimated, and the high-energy shoulder of the L 3 peak is not reproduced. The magnetic 
scattering is also well described, but it overestimates the intensity of the L2 peak. Both low- 
and high-energy shoulders in the L 3 peak are not present in the fit. We note that we could 
not obtain satisfactory simultaneous fits to the orbital and magnetic scattering using the 
other type of the orbital ordering {d 3x 2_ r 2 / d 3y 2_ r 2) . The fit corresponding to the crystal field 
described with X 400 = 5.6 eV and X 220 = -1.5 eV, reproduced the spectrum for the orbital 
scattering very well (having all the peaks as the fit on Fig. |2]with even slightly improved 
the L3/L2 ratio), but the magnetic scattering was poorly described, completely missing 
the first feature in the L 3 peak. We also note that the addition of a small orthorhombic 
(D 2 h) component, X 222 = -0.85 eV, to the crystal field related to the d x 2_ z 2 / d y 2_ z 2 type 
of the orbital ordering somewhat improved the L3/L2 ratio for the orbital scattering and 
broadened the previously underestimated width of the L3 peak for the magnetic scattering, 
but kept the rest of the spectra shown in Figs. 121 and El unchanged. 

In order to address the origin of the different spectral features, we have examined the 
effect of a reduction of the tetragonal part (X 220 ) of the crystal field, using different values 
of X 220 as shown in Fig. This approach is based on the fact that changes of the Jahn- 
Teller distortion yield corresponding changes to the tetragonal field. We conclude that the 
reduction of the Jahn- Teller effect results a reduced L3/L2 ratio and in the drastic reduction 
of the main feature at the L 3 edge. 

IV. TEMPERATURE DEPENDENCE 

The orbital order spectrum shown in Fig.Elwas obtained at a temperature of 63 K. We also 
made measurements of spectral features in the (^,|,0) orbital order reflection as a function of 
temperature, upon warming, to measure the change of the crystal field. These measurements 
were performed on station 5U1 and the energy resolution (AE) was deliberately reduced 
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FIG. 5: Variation of the Jahn- Teller effect through the reduction of the tetragonal crystal field: 
X 220 = 3 75 e y ( full black line ) 5 2 eV (dashed red line) and 0.1 eV (full green line). 

to approximately integrate the spectral feature over energy. These results are displayed in 
Fig. El 

Figure El shows that the scattered intensity at the L 3 edge is very sensitive to the Jahn- 
Teller effect whilst the intensity at the L2 edge is not. We define for the purposes of this 
discussion, the branching ratio (Fig. El middle panel) as the ratio of integrated intensities 
(over q) measured at energies of 643 eV and 653 eV for the L 3 and L 2 edges respectively, 
shown in Fig. El top panel. Therefore, this branching ratio gives an indication of the Jahn- 
Teller distortion. 

Upon cooling past T 00 the branching ratio increases, suggesting an increase in the Jahn- 
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FIG. 6: (a) Temperature dependence of the integrated intensity of the main features at the L3 edge 
(triangles) and at the L2 edge (inverted triangles), (b) The temperature dependence of the ratio 
of the integrated intensities L3/L2. (c) The temperature dependence of the integrated intensity of 
the (4,-4,5) magnetic reflection measured at 643 eV. 



Teller distortion. Approaching Tjy the integrated intensity for both spectral features rises 
sharply. However, the branching ratio also increases at this point indicating that the Jahn- 
Teller distortion becomes even more dominant. At Tn, as indicated by the measured mag- 
netic order parameter in the bottom panel of Fig. El both the integrated intensity and the 
branching ratio saturate. Upon further cooling only a slight decrease in the branching ratio 
is observed. 

At temperatures below Tjy, the orbital order correlations are approximately constant even 
as the anti-ferromagnetic spin correlations are gradually increasing. The dramatic increase 
in the intensity of the orbital ordering superlattice reflection at Tn is very unexpected. The 



stepwise increase in the integrated intensity suggests that the long-range orbital ordering 
interactions maximize at T/v- The temperature region between T/v and Too is one in which 
the orbital order correlations gradually increase with decreasing temperatures. These results 
demonstrate the strong interaction between the orbital and spin degrees of freedom. Is such 
behavior typical, or restricted to just Lao.sSri.sMnC^? Because our study is one of the 
first to directly measure both the orbital and magnetic correlations, and their temperature 
dependences, it is difficult to be certain. However, there are intriguing results obtained 
by other resonant X-ray studies at the K edge of other 3d transition metal compounds 
that suggest that such effects may be a general feature of strongly correlated systems. In 
LaMnOa RXS studies*^ at the Mn K edge observed scattering at the symmetry forbidden 
(300) position. Such measurements, indirectly attributed to orbital ordering via Jahn- Teller 
distortions, showed that the intensity increases dramatically at Too and thereafter remained 
almost constant until just above T/v where the integrated intensity increased by 50%. In 
KCUF3 a doubling of the integrated intensity of the orbital order reflection was observed 
upon cooling just above T/ v 36 i 37 . What is surprising about all of these results is that the 
increase occurs well below the orbital ordering transition temperature where the orbital 
order parameter is expected to be saturated. A number of different causes of this effect has 
been proposed. Ishihara and Maekawa^ suggested that such anomalous changes may signal 
a change in the type of orbital order. However, our energy resonance data do not show 
any dramatic changes in profile around T/v- Paolasini et a/.— suggested that the magnetic 
order is driven by the orbital order and that the spin exchange constants are determined 
by the relative orientation of the occupied orbitals. Alternatively, Binggeli and Altarelli^ 
have suggested that in KCUF3 the increase in the orbital order reflection at T/v may simply 
suggest a low temperature structural transition. In Lao.sSri.sMnC^ our results suggest that 
the intensity of the orbital order reflection at T/v is linked to an increase in the Jahn- Teller 
distortion in the absence of any structural transitions. 



V. CONCLUSIONS 



We have presented the energy dependence of the scattered intensity at a fixed wavevector 
of qoo =(i,i>0) and q AF ={\-\,\) through the Mn L 3 and L 2 edges. 

A good agreement was found between the experimental data and theoretical fits which 



included a realistic description of the spin direction. A single set of parameters describing 
the crystal field was used for both orbital and magnetic reflections. It is noted that there 
was a slight discrepancy in the branching ratio (L 3 /L 2 ) between data and theoretical results. 

The results of the fitting procedure demonstrate that the orbital ordering is predominantly 
of the d x 2_ z 2 1 'd y 2_ z 2 type. A satisfactory fit to the both data sets could not be obtained with 
the other proposed type of the orbital ordering (d 3x 2_ r 2 / d 3y 2_ r 2) . The inclusion of a small 
orthorhombic component of the crystal field, in the case of the d x 2_ z 2 / d y 2_ z 2 type of the 
orbital ordering, moderately improved the fits. Such results are consistent with the x-ray 
linear dichroism results of Huang et alM. In addition, the x-ray scattering results confirm 
the existence of long-range order. 

Further theoretical investigation of the orbital scattering indicates that the intensity of 
the L 3 edge is strongly sensitive to the Jahn- Teller distortion, whereas the L 2 edge is mostly 
sensitive to orbital order. We believe this is a general feature which can be used as a 
guideline in the interpretation of the orbital scattering in other manganite systems. We 
conclude therefore, in the specific case of Lao.sSr! 5 Mn0 4 presented here, that the system is 
highly influenced by Jahn- Teller distortions. 

Measurements of the temperature dependence of the ratio of the integrated intensities at 
the L 3 and L 2 edges show that there is a strong coupling between the Jahn- Teller distortion 
and the long-range magnetic order. These results are not consistent with a change of the 
orbital order type at T/v, as previously suggested^, but are suggestive of an increase of the 
Jahn- Teller distortions below T N . 
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